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Method and device for determining write parameters for recording information on a record 
carrier 



The present invention relates to a method of determining write parameters for 
recording information on a record carrier, said information being in the form of a multi- 
dimensional channel data stream to be recorded as a channel band of at least two symbol 
rows one-dimensionally evolving along a first direction and aligned with each other along a 

5 second direction. The present invention relates further to a method and a corresponding 
device for determining write parameters for recording information on a record carrier, said 
information being in the form of a channel data stream to be recorded as a channel band of at 
least one symbol row one-dimensionally evolving along a first direction, wherein the write 
parameters are determined by an iterative procedure. The present invention relates further to 

10 a recording method and a corresponding recording apparatus for recording information in the 
form of a channel data stream on a record carrier. Still further, the present invention relates to 
a computer program for implementing said methods and to a record carrier. 

The record carrier can generally be based on magnetic recording principles or 
on optica] recording principles. The further description focuses in more detail on an optical 

15 record carrier, which, however, does not exclude other types of record carrier. 

Generally, for the 2D optical recording channel as a whole, that is, being the 
combination of the write-channel at the transmitting end of the channel and the read-channel 
at the receiving end of the channel certain properties shall be achieved. One main goal is 
linearization of the channel. It is assumed that the read-channel is more or less fixed by the 

20 characteristics of central aperture (CA) detection, i.e. the detection mode commonly used in 
ID optical recording (see e.g. J. Braat, "Read-out of Optical Disks", in "Principles of Optical 
Disc Systems", Adam Hilger Ltd, 1985, pp. 7-87.). Non-linear characteristics of the read- 
channel have to be compensated by proper measures taken at the side of the write-channel: 
this is known as write-precompensation, implemented via a write-strategy. The channel 

25 symbols (bits or, more generally, M-aiy symbols) are processed through a so-called (non- 
linear) transmit filter to generate the parameters for the physical write-channel. In case of 
(small) deficiencies in the write-strategy leading to an incomplete linearization, the remaining 
non-linearities can be dealt with by a complementary receiving filter (via a non-linearity 
compensation). An appropriate write-strategy for 2D optical storage, in particular for ROM 
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media, is therefore desired. Further, also specific measures for a write-strategy in recordable 
and/or re-writable 2D optical storage are required. 

A write-strategy procedure that realizes a first desired property of the high- 
frequency (HF) signal values that are detected in 2D modulation on (quasi-) hexagonal two- 

5 dimensional lattices of bits has been described in European patent application EP 02 076 
255.5 (PHNL 020279). The "physical" detection is based on the principle of central aperture 
detection of the photon density incident on the photo-detector (PDIC). On the hexagonal 
lattice, a hexagonal cluster consisting of 7 bits, with one central bit and 6 (nearest) neighbour 
bits, is considered as a basic unit, also called symbol unit or bit cluster. The first desired 

1 0 property is that the HF signal values show a systematic roll-off with an increasing number of 
neighbour bits of the pit-type ("l"-bits): this property must hold for both possible bit-values 
for the central bit When this property is not satisfied (e.g. for the pit-bits), the problem of 
signal folding has to be dealt with, which implies that (part of) the HF signal values increase 
(instead of decrease) with an increasing number of neighbour pit-bits (when the central bit is 

15 of the pit-type); moreover, in the case of maximum signal folding, which occurs when the 
bit-cell for a pit-bit consists of 100% "pit-area", it implies that the HF signal for the all-land 
case is identical to the HF signal for the all-pit case (both behave as perfect mirrors). 

Signal folding typically occurs when the pit-bits are physically mastered (in a 
ROM disc) such that the pit-area covers a large fraction of or even the complete area of a bit- 

20 cell (which is a hexagon, the fundamental cell of the 2D hexagonal lattice). The elimination 
of signal folding was achieved through the writing of (relatively much) smaller pit-holes than 
the ones that are maximally possible: a quite convenient roll-off of the signal values is 
achieved for a duty factor of 50%, that is, the pit-hole covers about half of the area of the 
available hexagon. 

25 Apart from the above mentioned "first desired property", a second additional 

desired property to be realized through an extended write-strategy shall be achieved. A 
variety of "second desired properties" can be thought of. A very likely candidate is that the 
HF signal values exhibit a signal variation that is typical for a linear response. Many 
candidate bit-detection schemes expect a linear response; since this type of bit-detectors 

30 cannot deal with channel non-linearities, some kind of (possibly memory-less) non-linearity 
compensation (NLC) has to be included prior to equalization and bit-detection. There are two 
disadvantages in the use of such an NLC circuit: firstly, the (memory-less) NLC suffers from 
its limited accuracy; and secondly, assuming that the noise distributions are level- 
independent, it is advantageous in view of limiting the influence of noise to spread the HF- 
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signal values as much as possible over the available amplitude space: such a situation is not 
accomplished by the measures according to the "first desired property" since the signal levels 
for the "l"-bit are non-linearly compressed, resulting in non-equidistant signal levels prior to 
the NLC. The NLC operation will then result in noise distributions that are dependent on the 
5 signal level. Therefore, it is advantageous to incorporate a write-strategy that delivers "linear 
levels" (as linear as possible) at the output of the physical bit-detection on the photo-detector, 
prior to any signal processing: as a result, noise variances will be equal for each individual 
level. 

A so-called PIP TM (Pre-compensation Iteration Process) write-strategy for 

10 use in multi-level (ML) (one-dimensional) optical recording has been disclosed in WO 
01/57856, Therein, a dedicated write-strategy is based on a write-strategy matrix, which 
depends on the central symbol to be written, and a limited number of its neighbouring 
symbols. PIP is promoted as an adaptive ML write strategy, designed to remove the majority 
of non-linear channel effects. In particular, PIP makes data recovery more robust by reducing 

15 the overlap between the distributions of neighbouring signal levels, which is accomplished by 
decreasing the width of these distributions, and most importantly, by centering the 
distributions, making the levels of the multi-level system equidistant 

In a two-dimensional pattern, at a particular location exactly the same bit 
cluster (or symbol unit) and thus also the same cluster-class can appear. A write-strategy that 

20 is based on a write-strategy table or matrix, as disclosed in WO 01/57856 for a one- 
dimensional coding scheme, then yields exactly the same write-strategy parameter, for 
instance, exactly the same pit-hole radius. However, even if the bit cluster is identical, the 
bits surrounding the cluster can be different, so that individual bits of the bit cluster will have 
write parameters, e.g. pit-hole sizes, that are different from their nominal values. These local 

25 deviations from the nominal write parameters (pit-hole radii) will influence the optimal 

choice for the write parameters to be made at the central bit. This can be partly accounted for 
by extending the size of the write-strategy table, e.g. by inclusion of more rings or shells of 
neighbouring bits of the bit cluster. A full account of this "chain-effecf • of one bit 
influencing the choice of the write parameters of a neighbouring bit would be to have a very 

30 large write-strategy table, which is however impractical to work with. 

It is an object of the present invention to provide a a method and a 
corresponding device for determining the write parameters which can effectively be used for 
multi-dimensional coding schemes. It is a further object of the present invention to provide a 
method and a corresponding device for determining the write parameters which take into 
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account the above described "chain-effect", preferably avoiding the use of a very large write- 
strategy table or matrix. Furthermore, an appropriate recording method and recording 
apparatus, computer program and a record carrier using the invention, shall be provided. 

This object is achieved according to the present invention by a method as 
5 claimed in claim 1 wherein the write parameters for recording a pit-symbol of a symbol unit 
of said channel data stream, a symbol unit comprising a central symbol and a number of 
neighbouring symbols of which some are located on the same symbol row as the central 
symbol and others are located on neighbouring symbol rows, are determined under joint 
consideration of 

1 0 (i) the symbol value of the central symbol of the symbol unit; 

(ii) the symbol values of the neighbouring symbols of the symbol unit located in 
the same symbol row as the central symbol of the symbol unit; and 

(iii) the symbol values of neighbouring symbols of the symbol unit located in the 
symbol rows that are neighbouring the symbol row of the central symbol of the symbol unit. 

15 Contrary to the solution known from WO 01/57856 write parameters for 

recording a pit-symbol of a symbol unit depend not only on the neighbouring symbols in the 
same symbol row at which the symbol under consideration is located, but in addition, depend 
also on the neighbouring symbols in the symbol rows above or below the symbol row at 
which the symbol under consideration is located. Thus, symbol values of symbols in 

20 neighbouring symbol rows determine partly the write parameters of a symbol in a given row, 
in order to achieve characteristics of the HF-signal of said symbol in said given row. 

In an embodiment of the invention the write parameters are determined by use 
of a parameter table containing the write parameters for all possible classes of symbol units, 
from which the write parameters for recording a pit-symbol at the central symbol of the 

25 symbol unit are selected according to the actual said symbol unit For each value of the 

central symbol, and for each of the possible environments (of neighbouring symbols) of that 
central symbol, an entry in the parameter table (also called write-strategy matrix), which 
yields a set of (at least one) write strategy parameters) for the central symbol of the symbol 
unit under consideration to be used in the physical write-channel. Instead of a single write 

30 strategy matrix a set of write-strategy matrices (preferably in a rewritable system) can be 
used, for instance, when the write-strategy involves more than one physical parameter; 
another application with a set of write-strategy matrices relates to the case where each matrix 
is devised for one physical condition of the write-channel (for instance, like tilt of the writing 
laser-spot, relative to the disc). 
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5 

The object underlying the invention is further achieved according to the 
present invention by a method as claimed in claim 4 comprising the steps of: 

setting the write parameters for recording pit-symbols of said channel data 
stream to preliminary parameter values, 
5 - updating the preliminary parameter values by searching for the updated 

parameter values that best fulfil a predetermined criterion for the write parameters for 
recording of pit-symbols, said criterion being determined by the difference of HF-signal 
values, which will be determined by use of a channel model or obtained during read-out of 
pit-symbols recorded by use of the updated parameter values (that are updated in a previous 
1 0 iteration), and reference HF-signal values, 

iterating said updating until a predetermined condition is fulfilled. 
According to this second embodiment of the present invention a write- 
precompensation through an "on-the-fly" (iterative) computational procedure is proposed that 
operates sequentially for a sequence of channel symbols, preferably in (roughly) the order at 
1 5 which these symbols have to be written to the record carrier: the write parameters of a current 
channel symbol are derived from the (already determined) write parameters of (a limited set 
of) previous channel symbols together with the write parameters of (a limited set of) future 
channel symbols. For these future symbols, an average (preliminary) write parameter is set, at 
least in the first iteration of the described procedure. In next iterations, for the future channel 
20 symbols the write parameters that are obtained during the previous iteration can be used to 
update the current channel symbol. The write parameters for a cluster of symbols will thus be 
determined not only by the composition of that cluster, but also to some extent by the histoxy 
(memory) of the preceding sequence of channel symbols that leads to the considered cluster 
at a given position along the sequence of channel symbols, via the values of the write 
25 parameter that has been set at the channel symbols of that preceding sequence. 

The present invention is preferably applied for 2D optical recording where the 
information is in the form of a channel data stream to be recorded as a channel band of at 
least two symbol rows one-dimensionally evolving along a first direction and aligned with 
each other along a second direction, said two directions constituting a two-dimensional lattice 
30 of symbol positions. However, the invention is generally applicable, i.e. it can be applied also 
for multi-dimensional recording where data are arranged along a 3D (or theoretically higher- 
dimensional) array. 

Corresponding devices which are adapted for carrying out said methods are 
defined in claims 14 and 15. A recording method and a corresponding recording apparatus in 
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which pit-symbols are recorded by use of write parameters which are determined by a 
procedure as defined above are claimed in claims 16 and 17. A computer program 
comprising program code means for causing a computer to perform the steps of the above 
methods when said computer program is executed on a computer is defined in claim 18. 
5 Preferred embodiments of the invention are defined in the dependent claims. 

The predetermined criterion to be fulfilled for the write parameters is preferably determined 
by the sum of absolute values of the differences of the so-called "read-out" HF-signal values 
which are the HF-signal values obtained from or to be obtained from read-out, and the so- 
called reference HF-signal values or by the sum of squared differences of said read-out HF- 
10 signal values and said reference HF-signal values. Preferably, said sum comprises squared 
differences for all pit-symbols and non-pit symbols (or "land" symbols) in a particular 
symbol area and said sum shall be minimized during updating. 

According to further embodiments of the invention it is proposed that the 
predetermined condition is that the write-parameter for each pit-symbol has been updated for 
15 a predetermined number of times or that is has reached a value below a predetermined 

threshold value, so that the predetermined condition is a quality measure or figure-of-merit 
The reference HF-signal values are obtained from a hypothetical ideal signal 
which would result for a linear channel, that is a channel that can be represented by a linear 
(two-dimensional) impulse response. For the read-out HF-signals on the other hand, in 
20 practice a finite number of write parameters will be used for which the resulting HF-signal 
values are determined in advance based on a computational model that represents well the 
experimental signal generation in the read-channel. Then, with a proper minimization 
procedure as defined above, the set of write parameters that gives the best match between the 
"read-out" HF signals (obtained from the computational model for the signal generation) and 
25 the (linear) reference HF-signal values can be found. For a finite number of write parameters, 
such a minimization procedure could be solved with a dynamic programming approach just 
like the Viterbi algorithm as is used for bit-detection in the read-i ■ 
enormous complexity aspects related to an M-ary Viterbi in case of many possible pit-hole 
sizes for the write-parameter, it is preferred according to the present invention that a low- 
30 complexity and slightly sub-optimal optimization procedure is used for realizing the match 
that is concerned with the best set of write parameters that realize the closest match between 
the targeted linear HF-signal values and the "read-out" HF-signal values that can be either 
computed for the write parameter set derived from the computational channel model, or the 
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"read-out" HF-signal values that can be directly measured when the write-parameters have 
been adopted for writing the symbols in an iterative writing experiment. 

As already described above the "read -out" HF-signal values and the reference 
HF-signal values are determined on the basis of symbol units, also called bit or symbol 

5 clusters, each symbol unit comprising a central symbol and a number of neighbouring 
symbols, in particular a number of nearest neighbouring symbols surrounding the central 
symbol. Such a symbol unit can, for instance, be a hexagonal cluster comprising one central 
symbol and 6 surrounding symbols at a nearest neighbour distance. Alternatively, a squared 
cluster can be used comprising one central symbol and 4 nearest neighbouring symbols. For 

10 the case of the hexagonal cluster, two of the 6 nearest neighbour symbols are located in the 
same symbol row as the central symbol, and the other four nearest neighbour symbols are 
located in neighbouring symbol rows. 

Furthermore, for the iterative procedure of the second embodiment, it is 
preferred that the preliminary write parameter values for the pit-symbols set in the first step 

15 of the method are derived from a parameter table containing the write parameters for all 
possible classes of symbol units (with a central pit-symbol). Alternatively, assuming a 
channel with binary modulation, to all pit-symbols the same fixed write parameters could be 
assigned prior to the first iteration. 

According to a preferred embodiment of the second embodiment comprising 

20 the iterative procedure, the iterative optimization procedure according to the present 

invention is based on a sliding window approach according to which in said updating step of 
the iteration the write parameters of the pit-symbols to be updated are updated subsequently 
symbol column by symbol column for a number of symbol columns defining a detection 
window, wherein the detection window is shifted after each iteration by at least one column 

25 in the tangential direction of the broad spiral which comprises a number of bit-rows aligned 
with each other in the second direction, whereby the write parameters of symbols in a new 
column that enters the detection window after sliding are set to initial predetermined values, 
and wherein the iterations are repeated for a given column until said column is shifted outside 
of said detection window. This is a simple sequential procedure for updating the write 

30 parameters which can be easily implemented. 

The write parameters to be determined according to the present invention 
mainly depend on the type of record carrier to be used. For a read-only (ROM) record carrier, 
the pit-hole size needs to be determined which is realized during mastering by applying a 
certain laser intensity for illumination of a photo-resist layer. For a rewritable record carrier, 
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based on phase-change technology, a certain amorphous region is realized by a series of laser 
pulses at well defined laser powers. Thus, instead of pit-hole size the more direct physical 
parameters that yield a given pit-hole size can be determined, such as the characteristics of 
write pulses, in particular the number, the duration and/or the power level of a plurality of 

5 write pulses, or, in a more simple case, the power level of a single write pulse. 

A record carrier on which pit-symbols have been recorded by use of the 
method according to the present invention is defined in claim 19. It can be seen from the 
record carrier, for instance by use of a SEM, TEM or AFM image, whether the pit-hole sizes 
are all the same, independent of the bit cluster type or whether they are different depending 

10 on a cluster. In the latter case, it is even possible to distinguish between two cases: in a first 
case all clusters that occur for one cluster type lead to the same pit-hole size, which indicates 
that a write strategy matrix or table has been used. In a second case, clusters that occur for 
one cluster type may lead to slightly different pit-hole sizes because an updating strategy is 
used according to the present invention. In order to recognize whether the variation is random 

15 or according to a particular update strategy, the 2D correlation properties of the pit-hole sizes 
of a given cluster type can be evaluated as a function of its neighbouring symbols which then 
indicates that the second embodiment of the present invention has been used to determine the 
pit-hole sizes. 



20 

The present invention will now be explained in more detail with reference to 
the drawings in which 

Fig. 1 shows a block diagram of a general layout of a coding system, 

Fig. 2 shows a schematic diagram indicating a strip-based two-dimensional 
25 coding scheme, 

Fig. 3 shows a schematic signal-pattern for a two dimensional code on 
hexagonal lattices, 

Fig. 4 illustrates two types of bi-linear interferences in a hexagonal cluster, 
Fig. 5 shows a hexagonal bit cluster as used according to the present invention, 
30 Fig. 6 shows the HF-signal pattern as a function of the cluster type, 

Fig. 7 shows HF-signal patterns as a function of the cluster type for 2D 
modulation on a hexagonal lattice for various fixed pit-hole sizes, 

Fig. 8 shows a schematic diagram for an iterative method according to the 

invention, 
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Fig. 9 shows the basic cluster classes for a 7-bit hexagonal cluster. 
Fig. 10 illustrates the problem underlying the present invention, 
Fig. 1 1 illustrates the sliding window implementation of the present invention, 
Fig. 12 illustrates the method of the present invention in more detail, 
5 Fig. 13 shows a 7-bit hexagonal cluster with 12 surrounding bits beyond the 

first shell of nearest neighbour bits, 

Fig. 14 shows the cluster classes for the 7-bit cluster shown in Fig. 13, and 
Fig. IS shows a schematic diagram of another embodiment of method 

according to the invention. 

10 



Fig. 1 shows typical coding and signal processing elements of a data storage 
system. The cycle of user data from input DI to output DO can include interleaving 10, error- 
correction-code (BCC) and modulation encoding 20, 30, signal preprocessing 40, data storage 

IS on the recording medium 50, signal post-processing 60, binary detection 70, and decoding 
80, 90 of the modulation code, and of the interleaved ECC. The ECC encoder 20 adds 
redundancy to the data in order to provide protection against errors from various noise 
sources. The ECC-encoded data are then passed on to a modulation encoder 30 which adapts 
the data to the channel, i.e. it manipulates the data into a form less likely to be corrupted by 

20 channel errors and more easily detected at the channel output. The modulated data are then 
input to a recording device, e.g. a spatial light modulator or the like, and stored in the 
recording medium 50. On the retrieving side, the reading device (e.g. photo-detector device 
or charge-coupled device (CCD)) returns pseudo-analog data values which must be 
transformed back into digital data (one bit per pixel for binary modulation schemes). The first 

25 step in this process is a post-processing step 60, called equalization, which attempts to undo 
distortions created in the recording process, still in the pseudo-analog domain. Then the array 
of pseudo-analog values is converted to an array of binary digital data via a bit detector 70. 
The array of digital data is then passed first to the modulation decoder 80, which performs 
the inverse operation to modulation encoding, and then to an ECC decoder 90. 

30 In the European patent application EP 01 203 878.2 the 2D constrained coding 

on hexagonal lattices in terms of nearest-neighbour clusters of channel bits is described. 
Therein, it has been focussed mainly on the constraints with their advantages in terms of 
more robust transmission over the channel, but not on the actual construction of such 2D 
codes. The latter topic is addressed in the European patent application 02 076 665.5 (PHNL 
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020368), i.e. the implementation and construction of such a 2D code is described therein. By 
way of example, a certain 2D hexagonal code shall be illustrated in the following. However, 
it should be noted that the general idea of the invention and all measures can be applied 
generally to any 2D code, in particular any 2D hexagonal or square lattice code. Finally, the 

5 general idea can also be applied to multi-dimensional codes, possibly with isotropic 

constraints, characterized by a one-dimensional evolution of the code in a certain direction. 

As mentioned, in the following a 2D hexagonal code shall be considered. The 
bits on the 2D hexagonal lattice can be identified in terms of bit clusters. A hexagonal cluster 
consists of a bit at a central lattice site, surrounded by six nearest neighbours at the 

10 neighbouring lattice sites. The code evolves along a one-dimensional direction. A 2D strip 
consists of a number of ID rows, stacked upon each other in a second direction orthogonal to 
the first direction, and forming an entity over which the 2D code can evolve. The principle of 
strip-based 2D coding is shown in Fig. 2, Several strips that are coherently stacked one upon 
the other forms a broad two-dimensional band, which can be spiralled on an optical disc 

15 (such a band is also called a "broad-spiral")- Between successive revolutions of the broad 
spiral, or between neighbouring 2D bands a guard band of, for instance, one (empty) bit-row 
(filled with zero-bits, and are thus equal to land-marks) may be located. 

The signal-levels for 2D recording on hexagonal lattices are identified by a 
plot of amplitude values of the HF-signal for the complete set of all hexagonal clusters that 

20 are possible. Use is further made of the isotropic assumption, that is, the channel impulse 
response is assumed to be circularly symmetric. This implies that, in order to characterize a 
7-bit cluster, it only matters to identify the central bit, and the number of "1* -bits (or "0"- 
bits) among the nearest-neighbour bits (0, 1, 6 out of the 6 neighbours can be a "l"-bit). A 
"0"-bit is a land-bit in our notation. A typical "signal-pattern" is shown in Fig. 3. Assuming a 

25 broad spiral consisting of 1 1 parallel bit rows, with a guard band of 1 (empty) bit row 

between successive broad spirals, the situation of Fig. 3 corresponds to a density increase 
with a factor of 1.7 compared to traditional ID optical recording (as used in e.g. in the Blu- 
ray Disc (BD) format (using a blue laser diode with a wavelength of 405 nm, and a lens with 
a numerical aperture of NA=0.85). 

30 The basic origin of the channel non-linearity is the fact that the detected signal 

is related to the photon probability at the photo-detector. The photon probability is modeled 
(in scalar diffraction theory) as the squared modulus of the (complex- valued) photon wave 
function (which describes the interaction of the possibly aberrated wavefront of the photon 
with the phase- and amplitude-structures on the optical disc constituted by the pits and lands). 
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The relation between the photon wave function and the bits written on the disc is (at least) a 
linear one. Therefore, the relation between the photon probability function and the bits is (at 
least) a bi-linear one, the terminology bi-linear being used here to indicate a non-linearity of 
second order. 

5 For the sake of completeness, it is to be noted that the photon probability 

function is further integrated over the photo-detector: this yields the so-called central aperture 
signal, referring to the (mathematically equivalent) integration of the photon probability in 
the plane of the (exit) pupil. The channel model yields linear and bi-linear terms. Among the 
bi-linear terms, self-interference terms for each pit bit (close enough to the center of the 

10 illuminating spot), and cross-interference terms for each two-bit pair (with both pit-bits 

within the area of the illuminating spot) are obtained. These bi-linear terms are illustrated in 
Fig. 4. The cross-interference terms become quite small when the distance between both pit- 
bits of the pit-pair is larger than the nearest-neighbour distance of the hexagonal lattice 
(which is equal to the hexagonal lattice parameter denoted by a): it is therefore a good 

1 5 approximation (especially for intermediate densities) to consider only nearest-neighbour 
cross-interferences. 

If the interferences of the channel are further limited to the bits of the 7-bit 
hexagonal cluster, the HF signal can be modeled to a very good approximation as (assuming 
for simplicity a pit-depth with single-pass phase modulation equal to n/2, for maximum 

20 modulation in the central aperture signal; and assuming a fixed pit-hole radius for all pit- 
bits): 

HF = 1 - 4bo(lo - Soft) - 4n(l n - s,^) + 8nb 0 xo.n + 8p n x n ^ . 

This is essentially a 4-parameter model (one parameter for each term). The parameters and 

variables have the following interpretation: 
25 n: number of nearest-neighbours (of the central bit) being of the pit-type; 

b 0 : bit-value of the central bit ("1" for pit, "0" for land); 

lo: tap-value of linear interference for central bit; 

l n : tap- value of linear interference for (nearest) neighbour bit; 

So ; o: value for self-interference of central pit-bit; 
30 Sn^: value for self-interference of (nearest) neighbour pit-bit; 

Xo,n: value of cross-interference between central pit-bit and (nearest) neighbour pit-bit; 

x n ^: value of cross-interference between two (nearest) neighbour pit-bits (neigbors of the 

central bit), which are also nearest neighbours of each other; 

p n : number of (nearest) neighbour pit-pairs among the (nearest) neighbour bits. 
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The possible values of the parameter pn (and its average value<^>) are shown 
for different values of the number of nearest-neighbour pit-bits n in the following table : 



n 


# of Neighbour Pit-Pairs (pj 


<Pn> 


0 


0 


0 


1 


0 


0 


2 


0,1 


0.4 


3 


0,1,2 


12 


4 


2,3 


2.4 


5 


4 


4 


6 


6 


6 



5 The above equation holds only for a fixed pit-hole radius for all pit-bits. If a varying pit-hole 
radius is allowed, then a generalized form of the above equation should be used instead, 
which reads as: 

HF = l-4b 0 (lo[S 0 ]-MSo]) 

-4L 6 W bi(ln[Si]-MSi]) 

10 +8L 6 MbobiXo^[S 0 ,Si] 

+ 8 S 6 M bibi+iXn.nESi, Sm]) 

where the same terminology has been used as in the above equation, but with explicit 
reference to the pit-hole surfaces indicated by Si for the pit-surface of pit-bit i. The indexing 
system of the bits on the hexagonal cluster is shown in Fig, 5 (and it is assumed that bj is 

1 5 again identical to b i ). 

Fig. 6 shows the HF signal pattern for a = 165nm and pit-diameter b = 
122.5nm. From the plot, the 1 1 different signals according to the number of different p n 
parameters can be clearly observed. The average HF signal value (indicated by the solid line 
in Fig. 6) is obtained as the average over all clusters with a given value of n (between 0 and 

20 6). This average value is determined by the value of <Pn>, which is listed in the third column 
of the above table. Since x n ,n is a positive number, the graph shows an upward curvature for 
the higher values of n (for both cases: bo = 0 and bo = 1). Thus, in conclusion, there are two 
basic types of non-linearity in this model. Firstly, there is the non-linearity associated with 
the cross-interference Xn^j which is governed by p n . Secondly, there is the non-linearity 

25 associated with the cross-interference xo,n which depends on the number of pit-pairs that 

contain the central (pit-) bit and the pit-bits among the (nearest) neighbours of the central bit 
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(which number is defined as n): so the pre-factorof xo, n is proportional to the product nbo. 
Since xo,n is a positive number, the second type of non-linearity (the 4th term in the right- 
hand-side of the above equation) boils down to a different (less negative) slope of the linear 
interferences for the case where the central bit bo = 1 as compared to the case where bo 83 0. 

5 In the above mentioned European patent application 02 076 2SS.5 (PHNL 

020279) the use of a single radius for the pit-holes, irrespective of the type of cluster that the 
corresponding central pit-bit belongs to, is proposed as a satisfactory means against signal 
folding, in particular on a ROM disc. In Fig. 7 the HF-signal patterns for various (fixed) sizes 
(diameters b) of the mastered pit-holes can be seen for a hexagonal lattice parameter 

10 a=165nm and for a series of fixed pit-hole diameters of b=100nm, 120nm, 140nm and 

165nm. The HF-signals have been obtained through a scalar diffraction model tailored to the 
2D hexagonal lattice. 

Fig. 8 illustrates the basic principle of the method according to the second 
embodiment of the current invention using the iterative procedure. At the input, the 2D bit- 

1 5 pattern that has to be written to the disc is provided. For each bit location (denoted by 
coordinates (k, 1)), the information of the bit-cluster consisting of the central bit and its 
neighbour bits is retrieved. In a initialisation step the write parameters p°w of the (non-zero) 
bits bki are set to a preliminary value, e.g. a fixed write parameter (e.g. pit-hole size) or 
obtained from a table or matrix. Thereafter, these preliminary values are updated in an 

20 iterative procedure. 

The bit-cluster that is referred to may consist of the central bit plus a number 
of shells consisting of neighbour bits all at the same distance from the central bit. The 
simplest case is the one with only one shell (containing the nearest neighbour bits), yielding 
7-bit clusters. This one-shell case seems to be quite accurate for moderate to even relatively 

25 high recording densities in 2D optical recording. Therefore, it is treated in more detail in the 
following, as a representative but specific example. 

In principle, the write-strategy could be devised for any symmetry present in 
the read-out spot (like an elliptic shape). For simplicity's sake from now on isotropic (read-) 
channel characteristics, implying a read-channel with a circularly-symmetric symmetry, or at 

30 least a symmetry compatible with the hexagonal (rotation) symmetry of the 2D bit-lattice are 
considered. The basic (or independent) cluster classes for this case are now derived: a cluster 
class comprises all clusters that can be transformed one into another by means of rotation 
over 60, 120, 180, 240 or 300 degrees. It turns out that there are 28 of such independent 
cluster classes, 14 with the central bit value b 0 equal to 0, and 14 with b 0 equal to 1 
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(considering only a non-zero pit-hole radius for pit-bits, having bit-value bo m 1). These basic 
cluster classes are denoted in Fig. 9 as PAT-01, PAT-02, PAT-14. In order to describe the 
different cluster classes, we have adopted the convention as shown in Fig. 9. For each cluster 
class, its multiplicity (denoted by xi) is indicated by the number "i" which is the number of 

5 clusters that belongs to a given cluster class. It should be noted that the (rotation variants of) 
classes PAT-08 and PAT-09 can be transformed one into another by point inversion (with the 
centre of inversion located in the centre of the cluster). So, if the inversion symmetry is added 
then the number of distinct cluster classes reduces to 13 (PAT-08 and PAT-09 becoming 
degenerate). A next reduction in number of distinct classes is possible if only next-neighbour 

10 non-linearities for xo,n are taken into account. Then, classes PAT-03 and PAT-04 become 
degenerate; the same holds for classes PAT-10 and PAT-1 1. Thus the number of distinct 
classes has become equal to 1 1 . A further (and still more severe) reduction to only 7 distinct 
classes is possible if only the number of neighbour pit-bits n as a relevant parameter is 
considered. 

1 5 The problem underlying the present invention shall now be illustrated with 

reference to Fig. 10. At location (k,l) for the two situations indicated by a circle exactly the 
same cluster CI, C2 and thus also the same cluster-class is found as shown in Fig. 10a and 
Fig. 10b. A write-strategy that would be based on a write-strategy table or matrix would yield 
exactly the same write parameters. 

20 However, in situation (2) shown in Fig. 10b, the pit-bit b 5 2 at 7:30 in the circle 

is surrounded by three + one (the central bit of the cluster C2) pit-bits, whereas in situation 
(1) shown in Fig. 10a, the same pit-bit bsi is only surrounded by one+one pit-bit. A situation 
with no pit-neighbours for that bit outside of the circle determining the cluster-class is also 
possible. 

25 The bit bsi and b 5 2, respectively, will thus have a different pit-hole size, or 

more generally, different write parameters for both situations. These different sizes will 
influence the optimal choice for the pit-hole size to be made at the central bit bw. This can be 
partly accounted for by extending the size of the write-strategy table, e.g. by inclusion of 
more rings or shells of neighbouring bits (making the circle drawn larger and larger). A full 

30 account of this "chain-effect" of one pit-hole influencing the choice of a neighbouring pit- 
hole would be to have a very large write-strategy table, which is however impractical to work 
with. 

It is thus proposed according to the present solution to perform an "on-the-fly" 
optimization of the pit-hole sizes, taking into account the above "chain-effect". Instead of pit- 
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hole sizes (for ROM), it is also possible to optimize any set of parameters on which the write- 
channel (e.g. a set of laser-pulses with for each pulse a certain duration and a certain laser 
power, for phase-change recording) may be based, such as the power level or number of 
write pulses. 

5 For the following explanation it shall be supposed that there are L possible 

values for the pit-hole size of a pit-bit, that the memory of the read-channel (that is, the extent 
of the ISI) amounts to M fish-bones (one column or zig-zag pattern of channel symbols in the 
radial direction in the channel strip shown in Fig. 1 1) at each side of a current fish-bone (total 
memory is 2M) and that there are N row bit-rows within one channel strip. Then, the 

10 optimization of the write-strategy proposed according to the present invention, based on a 
quantitative figure-of-merit, is a dynamic programming problem, just like the standard 
Viterbi-problem. Optimization then means finding the best path (with the lowest cost, or 
lowest value of the Figure-of-Merit) through a trellis of states: the number of different states 
amounts to IA2 M Nrow). This number is a maximum, since some states may be forbidden, 

1 S because land-bits have always a zero pit-hole size which does not need optimization. 

However, instead of "exactly" solving the above optimization procedure with 
the use of the Viterbi algorithm, it is proposed to use an iterative procedure. As an example 
the optimization parameters are the pit-hole sizes of all pit-bits. The optimization criterion 
(or, Figure-of-Merit, FoM) reflects the intention to adopt the overall response of the channel 

20 (write-channel and read-channel) to a certain specified target-response. A convenient target- 
response (in view of a bit-detector) could be a linear one. A preferred embodiment of a FoM 
to be used is: 

FoM = Skj [HFchannd (bki + neighbouring bits) - HF target (bw + neighbouring bits)] 2 . 

25 

In the above equation, the first set of HF-values are the so-called "read-out 9 * 
HF-signal values, and the second set of HF-values are the so-called "reference" HF-signal 
values. 

The figure-of-merit should be made small enough by the write-strategy 
30 optimization. The FoM is the sum of squared values of the deviations of a target signal 

waveform HF<target> (which can be a linear target, but it can also be a non-linear target in 
order to use a larger portion of the signal amplitude range in the area of signal overlap in the 
2D signal pattern, said overlap occurring between the signal levels for the clusters where the 
central bit is a "0" and the signal levels for the clusters where the central bit is a "1") 
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subtracted from the signal waveform that results from a computational model (or, 
equivalent^, could be measured experimentally) for a given set of pit-hole sizes (that have 
been defined for the central pit-bit and its neighbour pit-bits). 

Optimization is done on the basis of a sliding window as shown in Fig. 1 L A 
5 rectangular window W is chosen as a practical example. The window W comprises all bit- 
rows in the two-dimensional broad spiral. Its lateral (or tangential) extent is a number of 
(N+l) fish-bones. 

Supposing the (one-sided) tangential extent of the ISI (in one direction) of the 
read-channel amounts to M fish-bones. Then, as shown in Fig. 1 la for time moment "k", for 

10 the pit-bits in the M fish-bones to the right hand side of the window W are set to an initial 
value; this initial value can be a constant value, or a value derived from some write-strategy 
table. First, the most right-positioned fish-bone F° in the window W at moment "k" gets 
updated values for its pit-hole sizes (denoted by the array S 0 ^ hereby using the pit-hole radii 
of M fish-bones at its right hand side, and the pit-hole radii of M fish-bones at its left hand 

15 side. Then the same updating procedure is used for the 2nd most right-positioned fish-bone 
F 1 in the window at moment "k"; and so-forth, until the left-boundaiy of the window W has 
been reached. For a general fish-bone inside the window W, the pit-hole radii of pit-bits to its 
right hand side have been updated before (for the same position of the window, but during 
optimization of the pit-hole sizes in a preceding fish-bone), whereas the pit-hole radii of pit- 

20 bits to its left hand side have been updated for the previous position of the window W (at 
moment "k-1" if the current moment is "k"). So, for a given position "k" of the window W, 
the pit-hole sizes of the fish-bones denoted by the arrays S° k , S l k -i up to S N k . N are 
successively being determined through this optimization procedure that proceeds fish-bone- 
by-fish-bone. 

25 After completion for the current position, the window W shifts one fish-bone 

to the right, and the optimization procedure starts all over again. This situation is shown in 
Fig. 1 lb for time-moment "k+1". 

The pit-hole sizes of a given fish-bone are thus iteratively updated according 
to the above procedure, with the number of iterations (or updates for a given fish-bone) equal 

30 to N+l . Only pit-hole sizes for non-zero bits are updated, the "0"-bits remain equal to the 
surrounding land. In Fig. 12 three situations for different pit-hole sizes to be optimized are 
shown. The central pit bo is the one of which the pit-hole size has to be updated. The pits b u 
have already been updated for the current position of the window; the pits b n have not yet 
been updated for the current position of the window, but have been updated at a previous 
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position of the window, or for the first iteration of the procedure, they have a pit-hole size 
from an educated guess, for instance from a write-strategy table. So, the pit-hole sizes from 
all the pits in a cluster of bits centered around the pit-bit bo to be updated are known. A 7-bit 
cluster has been used for convenience, but it might be any other (larger) cluster. 

5 The pit-hole size of the central pit-bit is also known, e.g. its value from the 

previous position of the window, from the educated guess or from a table. With all this 
knowledge of pit-hole sizes, the pit-hole size of the central pit-bit bo can be updated. For 
instance, 2N P +1 values of its pit-hole size can be considered, with a resolution or step size 
equal to "delta", centered around the previous value of the pit-hole size. For each of these 

10 candidate pit-hole sizes, or for a limited subset of these candidate pit-hole sizes, centered 
around the previous value of the write parameter, the FoM, in fact the terms that depend on 
the value of the pit-hole size of the pit-bit that is varied, are evaluated. For the 7-bit cluster, 
these are the 7 HF-signal values around and including the central pit-bit bo. 

The (most probably linear) target that has been set (i.e. HF<target>) is known. 

15 The actual HF-signal value HF<channel> for each of the 7 locations around and including the 
central pit-bit are then derived by a channel model that depends on linear and non-linear ISI 
coefficients that are explicitly dependent on the pit-hole sizes of the pit-holes at the pit-bits of 
interest 

In the following first embodiment of the present invention using a parameter 
20 table (write strategy matrix) for the determination of the write parameters shall be explained. 
The case with 14 distinct cluster classes will be considered as a practical example for the 
subsequent description. Further, a ROM optical disc for 2D optical storage will be 
considered, where the write strategy matrix contains as a write-strategy parameter for each 
class the area of the surface of the corresponding pit-hole (if the central bit is a pit-bit). The 
25 write-strategy matrix is derived in the following optimization procedure. For a given write 
strategy matrix, a figure-of-merit (FoM) is defined based on the squared value of the 
difference between the read-out HF signal and the desired HF signal (which is due to the 
target linear interferences only); this squared difference is derived for each cluster class 
(averaged over all possible values for the 12 surrounding bits b 7 , b 8 , bis which determine 
30 the actual surfaces to be used for the neighbour pit-bits in the cluster class) and the sum over 
all classes multiplied with the corresponding multiplicity factor yields the final value for 
FoM, given by: 

FoM = X £m, X £(^(& 0 ;c/a*O 

/=1 bo~0 67=0*8=0 418=0 
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where HF(bo; class*) has implicit dependence on the 12 next-neighbouring bits Q>% bs, bis) 
and where the desired (target) linear HF signal is given by: 
HFiin (bo; classi) - 1- bo Co - nj C|. 



5 The coefficients Co and ci are the central tap and neighbour tap coefficients of 

the desired linear response. The parameter nj is the number of neighbour pit-bits for class i. 
The figure-of-merit is a statistical average computed in a deterministic way by averaging 
over all possible cluster classes (2 7 distinct cases) and all possibilities of their surrounding 
bits (2 12 distinct cases). Via the introduction of the cluster classes, the number of 
10 computations is largely reduced, however, the resulting figure-of-merit remains exactly 
identical. 

Fig. 13 shows the 7 bits of a cluster together with the 12 surrounding bits l>7, 
bs, bis. Fig. 14 illustrates for a given cluster how the cluster classes of each pit-bit is 
determined (from the bit-values of its neighbour bits): for the central bit b 0 of the cluster, the 

15 bit-values of the cluster-bits are sufficient; for each of the neighbour bits in the cluster (bits 
b2, b3, b 6 ), 3 surrounding bits are required in order to uniquely determine its cluster class. 

Given a performance criterion like the above figure-of-merit, the search for the 
optimum write strategy matrix is simply an optimization procedure in a Nd-dimensional 
space (with Nd the number of cluster classes used). Brute-force search procedures are not 

20 suitable because of the large dimensionality of the problem. Any sub-optimal optimization 
procedure (like steepest-descent etc.) will be sufficient. 

A practical optimization procedure, that has been applied, considers a fixed 
number of distinct pit-hole surfaces. For each surface Si, the coefficients for linear 
interferences Io[S;] and l n [Sj] are computed; also all coefficients for self-interferences So,o[S;] 

25 and s Iwl [Si], and all combinations for pit-hole surfaces Si and Sj of the cross-interferences 
xo,n[Si; Sj] and x^pSi; Sj] have to be available too. From these parameters, the HF signal for 
any cluster with any possible pit-hole sizes out of the set of available pit-hole sizes can be 
computed. For each cluster, it is evaluated whether it is beneficial to decrease or increase the 
pit-hole-size of the central pit-symbol by a small amount (the step size used in the 

30 optimization); this procedure is performed for all clusters, after which it can be repeated in a 
number of subsequent iterations of the optmization procedure. 

The basic principle of this embodiment is illustrated Fig. 15. At the input, the 
2D bit-pattern that has to be written to the disc is provided. For each bit location (k,l), the 
information of the bit-cluster consisting of the central bit and its neighbour bits is retrieved. 
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Next, it is analysed to which cluster class (denoted by pi) the current cluster belongs. For the 
identified cluster class (at location (k,l)), the corresponding write parameters, for instance the 
size of the pit-hole for ROM mastering, is obtained from the write strategy matrix S. For a 
ROM system this matrix S contains the pit-hole sizes for the various basic cluster classes. 
5 This procedure is carried out for all bit-locations (k,l) when occupied by a pit-bit. It should be 
noted that it has been assumed here that a land-bit remains virginal, i.e. no pit-hole at all is 
mastered. 

In a practical implementation a lattice parameter a=165nm for NA=0.85 and 
X=405nm is considered, yielding a capacity increase over the BD-format of 1.4x. A write- 

10 strategy matrix has been derived in a simulation set-up based on scalar diffraction 

computations. In the optimization procedure, 40 possible pit-hole sizes with a surface ranging 
equidistantly from 0 up to 0.25rca 2 has been allowed. It has been observed that the linearized 
levels are very close to the target levels, illustrating the adequate performance of the write- 
strategy in linearizing the read-out signals. 

15 With a larger modulation (minimum modulation level at 5%), larger pit-hole 

sizes are required on the average than for the case of a smaller modulation. It is further noted 
that the average pit-hole size is markedly smaller than in the original situation without write- 
strategy (with the fixed pit-hole diameter of 122.5nm). The average pit-hole diameters for the 
two cases are: 97.8nm and 106.0nm. At this resolution in pit-hole radii (40 equidistant steps 

20 from zero to maximum pit-hole surface), cluster classes with the same value for n (the 

number of nearest neighbour pit-bits) show identical pit-hole radii, which makes the number 
of distinct entries in the write-strategy matrix even lower (reduction from 14 down to 7). 

In a further implementation a lattice parameter a=138nm for NA=0.85 and 
X=405nm is considered, yielding a capacity increase over the BD-format of 2x. The same 

25 procedure as described in the above paragraph is repeated, i.e. with pit-hole diameter 

b=102.5nm for the two cases with 15% and 5% minimum modulation level. The average pit- 
hole diameters for the two cases are: 83.6nm and 90.0nm. Similar observations have been 
made as in the previous section, apart from two aspects: a) the pit-radii do not "cluster" 
according to the number of nearest pit-neighbours n; and b) the maximum pit-hole size scales 

30 as the size of the hexagonal bit-cell (with a factor of about 1 .41) and the minimum pit-hole 
size scales down less fast (with a factor of 1.33). 

For the ease of the description, the interferences have been limited to the first 
shell of nearest neighbours. To obtain a well-suited write-strategy for an increased capacity 
(like 2x BD), however, it is preferred to include at least the 2nd shell, which will lead to a 
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larger number of entries in the write-strategy matrix. The average pit-hole size is markedly 
smaller than in the original situation without write-strategy (with fixed diameter b=102.5nm, 
which may be favourable in view of the proximity effect in electron beam recording (EBR). 
The invention is not limited to the 2D hexagonal lattice, but can also be 
5 applied to any type of 2D bit-lattice. Further, the invention is not limited to write-strategies 
that account for the interferences from the nearest neighbours (or first ring (or shell) of 
surrounding bits), but can be generalized to other (larger) sets of neighbouring bits. In an 
electron beam recorder, a write-strategy for minimization of the proximity effect (long-tail of 
write-impulse due to back-scattered electrons, can be up to l^m) may be required (certainly 

10 at high densities). The present proposal for linearization of the 2D channel may lead to a joint 
write-strategy, satisfying both purposes (that is linearization of the overall channel, and 
reduction of proximity-effects). 

The "desired property" to be realized by the invention can be different from 
"linearization". A possible candidate is to achieve a situation where the read-levels show a 

15 (much) reduced dependence on the two bits of the bottom row of a 7-bit hexagonal cluster. 
Such a situation is largely advantageous for a stripe-wise Viterbi bit-detector as described in 
European patent application 02292937 (PHNL021237EPP), in which a 2-row Viterbi- 
detector is used which shifts row-by-row from top to bottom of the broad spiral (consisting of 
more than 2 rows). In a known bit-detector the detector produces hard-decision or soft- 

20 decision information, and is used in an iterative way: this iterative processing is needed 

because the detector does not know in the first iteration the (probability of the) bits in the bit- 
row below a given (current) position of the stripe of the 2-row Viterbi-detector. Via a proper 
write-strategy, it is intended to reduce the impact of these two bits in the row "below": in this 
way, the number of iterations required in the stripe-wise Viterbi detector will be reduced. As 

25 such, this procedure is then a substitute procedure for a two-dimensional version of decision- 
feedback equalization at the side of the read-channel. Further, a combined solution with a 
transmit filter and a receive filter is possible. 

It is also possible to limit the resolution in pit-surface to only a few (e.g. 3) pit- 
sizes. Another viewpoint is to make the resolution in hardware not worse than the resolution 

30 that can be handled in the write process (e.g. the statistical variance in pit-size induced by a 
laser-beam recorder (LBR) or electron-beam recorder (EBR) ). 

According to the present invention write parameters for recording a pit-symbol 
of a symbol unit depend not only on the neighbouring symbols in the same symbol row at 
which the symbol under consideration is located, but in addition, depend also on the 
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neighbouring symbols in the symbol rows above or below the symbol row at which the 
symbol under consideration is located. Thus, symbol values of symbols in neighbouring 
symbol rows determine partly the write parameters of a symbol in a given row, in order to 
achieve characteristics of the HF-signal of said symbol in said given row. 
5 According to preferred embodiment of the present invention a solution is 

proposed to perform an "on-the-fly" optimization of the write parameters, in particular pit- 
hole sizes, for recording pits on a record carrier, taking into account the above described 
"chain-effect" where the size of one pit-hole at a given pit-bit is influenced by the chosen 
sizes of many neighbouring pit-holes. Instead of pit-hole sizes (for ROM), any set of 
10 parameters on which the write-channel (e.g. a set of laser-pulses for phase-change recording) 
may be based, can be optimized. 



